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Abstract 

In this article, we establish a probabilistic representation for the second- 
order moment of the solution of stochastic heat equation in [0, 1] xR'', with 
multiplicative noise, which is fractional in time and colored in space. This 
representation is similar to the one given in [8] in the case of an s.p.d.e. 
driven by a Gaussian noise, which is white in time. Unlike the formula 
of [3, which is based on the usual Poisson process, our representation is 
based on the planar Poisson process, due to the fractional component of 
the noise. 
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1 Introduction 

The classical Feynman-Kac (F-K) formula gives a stochastic representation for 
the solution of the heat equation with potential, as an exponential moment of a 
functional of Brownian paths (see e.g. [S]). This representation is a useful tool 
in stochastic analysis, in particular for the study stochastic partial differential 
equations (s.p.d.e. 's). We mention briefly several examples in this direction. 
The F-K formula lied at the origin of the existence, large time asymptotic and 
intermittency results of [4J and 5 , for the solution of the heat equation with 
random potential on Z'^, respectively R''. The same method, based on a dis- 
cretized F-K formula, was used in [IT] for obtaining an upper bound for the 
exponential behavior of the solution of the heat equation with random potential 
on a smooth compact manifold. The technique of [J was further refined in [6] 
in the case of parabolic equations with Levy noise, for proving the exponential 
growth of the solution. The F-K formula was used in 9J for solving stochastic 
parabolic equations, in the context of white noise analysis. A F-K formula for 
the solution of the stochastic KPP equation is used in [13] for examining the 
asymptotic behavior of the solution. 
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The present work has been motivated by the recent article [5j , in which the 
authors obtained an alternative probabilistic representation for the solution of 
a deterministic p.d.e., as well as a representation for the moments of the (mild- 
sense) solution of a s. p.d.e. perturbed by a Gaussian noise F, with "formal" 
covariance: 

More precisely, in {F{h),h G V} is a zero-mean Gaussian process with 
covariance E{F{h)F{g)) — {h,g)-p, where V is the completion of {l[o,t]xA;i G 
[0, 1], A 6 BbQS.'^)} with respect to the inner product (•, •)-p given by 

(l[o.t]xA,l[o.s]xs)p = As) / / f{x-y)dydx. 

J A Jb 

(Here, BbiW^) denotes the class of bounded Borel sets in M''.) 
In the particular case of the stochastic heat equation: 

du 1 

— = -Au + uF, t>0,xeR'^ (1) 
the representation for the second-moments of the (mild-sense) solution u is: 



w{t - rjv, , Bl^^ )wit ~ , Bl^^ ) n f{Bl^ - B^ ) 

i=i 

^ (2) 

(with the convention that on {Nt — 0}, the product is defined to be 1), where 
B^ = {Bl)t>Q and B^ ~ {B^)t>o are independent d-dimensional Brownian 
motions starting from x, respectively y, N = {Nt)t>o is an independent Poisson 
process with rate 1 and points ti < T2 < . . and 

■w(t,x)^[ pt{x -y)uo(y)dy, where pt(a;) —-^-T7^e"l"=l'/(^*). 

Note that the representation ([2]) docs not rely on the entire Brownian path, 
but only on its values at the (random) points ri , r2 , . . . , tj^^ . This property has 
allowed the authors of ^ to generalize the representation to a large class of 
s.p.d.e.'s, including the wave equation. 

To see where the idea for this representation comes from, we recall briefly 
the salient points leading to ([2]). If the solution of ([T]) exists, then it is unique 
and admits the Wiener chaos expansion: (see e.g. Proposition 4.1 of [8]) 

oo 

Ut,x = W{t,x) +^In{fn{-,t,x)), (3) 

71=1 

where /„(/„(•, i,a;)) J^^^^^^^^aj^ fniti,xi, . . . ,t„, Xn)dFt^_^^ . . . dFt^^^^, and 
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fn G V®" is a symmetric function given by: 

1 " 



(4) 

Here, p is a permutation of {l,...,7i} such that tp(i) < tp{2) < •■• < tp{n), 
tp{n+i) = t and a;p(„_|_i) = x. By the orthogonahty of the terms in the series ([3]), 

oc 

E[ut,xUt.y] = x)w{t, y) + ^ a;, y), 



n=l 



where 

Jn{t,x,y) 

F{ti, . . . , t„) 



«!(/„(•, i, x), /„(•, i, = / F{ti, . . . ,tn)dt, and 



/(xj - yj)dydxdt. 

Here, we denote t = {h, . . . ,t„),x = {xi, . . . ,a;„), y = (yi, ...,?;„) and T„(i) 
{(ti,...,t„);0<ii < ... <t„ <t}. 

A crucial observation of [8 is that, for any F : [0,t]" M+ measurable. 



/ F{ti, . . . , t„)dt = e*E^[F{t - r„, . . . , i - ti)1{^^^„}] 

JT„(t) 



(5) 



This is a key idea, which yields a probabilistic representation for Jn{t,x,y), 
based on the jump times of the Poisson process. This idea is new in the lit- 
erature, although it appeared implicitly in the earlier works [10], [13], [T6] . 
Secondly, for any < ii < . . . < i„ < i, F{t — tn, ■ ■ ■ ,t — ti) is represented as: 



F{t-tn,...,t-ti) = 



E 



,{t - t„, Bl)w{t - i„, Bl) Yl f{Bl^ - Bl) 



(6) 



Relation ^ follows from these observations, using the independence between 
N and 51,52. 

In this article, we generalize these ideas to the case of the stochastic heat 
equation driven by a fractional-colored noise. More precisely, we consider the 
following equation: 



-Au-fMoM/, t e [0, l],a; e M'' 



du 

at ^ 2 



(7) 
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where uq G Cb(R'*) is non-random, o denotes the Wick product, and is a 
Gaussian noise whose covariance is formally given by: 

with 



Here, H G (1/2, 1), an = H{2H — 1) and /i is a tempered measure on R''. (Note 
that r]{t, s) is the covariance kernel of the fractional Brownian motion.) 

The noise W is defined rigourously as in ,2^, by considering a zero- mean 
Gaussian process W — {W{h); h G HV} with covariance 



E{W{h)W{g))^{h,g) 



nv, 



where HV is the the completion of {l[o,t]xA; ^ G [0, 1], A G Bb(M.'^)} with respect 
to the inner product (•, ■)n'P given by: 

{^[o,t]xAA[a.s]xB)nv = / / / viu,v)f{x - y)dydxdvdu. 
Jo Jo JaJb 

(See also [7] for a martingale treatment of the case 'ri{t,s) = 5o{t — s), which 
corresponds to H = 1/2.) 

As in ^3J , the solution of equation ([7]) is interpreted in the mild sense, using 
the Skorohod integral with respect to W. More precisely, an adapted square- 
integrable process u ~ {ut^x', {t,x) G [0, 1] x R'*} is a solution to ([7]) if for any 
{t,x) G [0,1] X M'', the process {ptsix - y)us,yl[o^t]{s); (s,y) G [0,1] x M'*} is 
Skorohod integrable, and 

ut,x = PtUo{x) + / pt-s{x - y)us.yl[o,t]{s)SWs,y. 
Jo Jm-i 

The existence of a solution u (in the space of square-integrable processes) 
depends on the roughness of the noise, introduced by H and the kernel /. 
We mention briefly several cases which have been studied in the literature. If 
/ = (5o, the solution exists for d ^ 1,2 (see [12]). If fix) = Oti \xi\^"'-^, 
the solution exists for d < 2/{2H — 1) + J2i=i Hi (see [H]). Recently, it was 
shown in [3] that if / is the Riesz kernel of order a, or the Bessel kernel of order 
a < d, then the solution exists for d <2 + a, whereas if / is the heat kernel or 
the Poisson kernel, then the solution exists for any d. 

Our main result establishes a probabilistic representation for the second 
moment of the solution, similar to and based on the planar Poisson process. 

Theorem 1.1 Suppose that equation ^ has a solution u = {ut.x'i (t, x) G 
[0, l]xM''} in[0,l]xR'^. Then, for any {t,x) G [0, l]xM'' and{s,y) G [0,l]xR'', 

E[ut^xUs,y] = w{t, x)w{s, y) + e*^ 
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3 = 1 

where the sum is taken over all distinct indices ii, . . . , z„, 

• B^ = (i?t^)fG[o.i] and B^ = (Bf^)te[o.i] are independent d- dimensional stan- 
dard Brownian motions, starting from x, respectively y, 

• N = {Nt^s){t.s)e[o,i]^ is an independent planar Poisson process with rate 1 
and points {Pi}i>i, with Pi = {Ti,pi), 

• ^ii....,i„ (i, s) is the event that N has points Pi-^,...,Pi^ in [0,t] x [0, s], 
r* = max{rii , . . . , } and p* = max{pjj ,pi^}. 



21 ,...,in 



2 Proof of Theorem 11.11 

We begin by recalling some basic facts about the planar Poisson process. If 
N = {Nt.s)(t.s)e[o,i]^ is a 2-parameter process, and R — (a, b] x (c, d] is a rectangle 
in [0, we define Nr = Na,c + Nf,^d ~ Na,d - N^^c- 

We say that N — (-/Vt,s)(t,s)6[o,i]2 is a planar Poisson process of rate 
A > 0, if it satisfies the following conditions: 

(i) N vanishes on the axes, i.e. Ntfi = iVg.t = for all t G [0, 1]. 

(ii) Nfi has a Poisson distribution with mean A|i?|, for any rectangle i?; 
(Hi) Nji-^, . . . , Nr^. are independent, for any disjoint rectangles Ri, . . . , Rk- 

The following construction of the planar Poisson process is well-known (see 
G-g- llj)- Lst X a Poisson random variable with mean A, and {Pi}i>i be an 
independent sequence of i.i.d. random vectors, uniformly distributed on [0, 1]^. 
We denote Pi = [Ti.pi). For any {t,s) G [0,1]^, define 

X 

Nus^^I{T^<t,p^<s}. 
i=l 

Then N = (-/Vt,s)(t.s)G[o.i]2 is a planar Poisson process with rate A. 

For any n > 1 and for any distinct positive integers ii, . . . , i„, let Aii,...,i„ (i, s) 
be the event that N has points Pi^ , ■ ■ ■ , Pi^ in [0, t] x [0, s] . Then 

{Nt,s=n}^ y (8) 

2i,...,i7i>l distinct 

The following result is probably well-known. We include it for the sake of 
completeness. 

Lemma 2.1 The conditional distribution of 

{{t-Ti^,s- piJ,...,{t-Ti^,s- p^J) given Ai^^,,,j^{t, s) 
is uniform over {[0,t] x [0, s])". 
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Proof: Let / = {ji, . . . ,i„} and i* = maxjii, . . . ,z„}. Then 

Ai{t,s) = 



fe>n,/C{l,..., 

and 



U {X = fc}f| (f|{T, <t,p, <4] fl [ fl {r, Xorp, >s} 
C{l,...,/c} Vie/ / yi<fc,i^/ J 



A:! 

fe>nVi* 



For any Borel sets Fi, . . . , r„ in [0, t] x [0, s], we have 

ij = l / fc>nVi* ■ \j = l 



and hence 



\ 1 " 



pyf]{it-n^,s-p,^)er,}\Ai{t,s) 
□ 

As a consequence, for any measurable function F : ([0,t] x [0, s])" R-|_, 
£;[F(i- Til, s-pii,...,t s-Pj,J I s)] = (9) 

-A--/ / F(ii,si,...,i„,s„)dsdt, 
[tsr J[o,t]" ^[0,^]" 

where t = (ii, . . . , i„) and s — (si, . . . , s„). 

Suppose that A = 1. Then (ts)" = n! e*''P(iVt,3 = n) and 



■ / / -^"(^1,31, . . . ,t„,s„)dsdt. 



E[F{t - 3 - Pi^, . . . ,t - Ti^^, S ~ p^J I Aii,...^i„(t, s)] = 

1 

nle^'^PiNt^s = n) J[„ f]„ J[o ,,y. 

Using ([8]) and ([9]), we obtain that for any F : {[0,t] x [0, s])" — > R+ measur- 
able, 

/ / F{ti,si,...,tn,Sn)dsdt= (10) 

J[0,t]" J[0,s]" 
ii^...,in>l distinct 

Relation (fTO|) is the analogue of ((S]), needed in the fractional case. 

Suppose now that equation ([7]) has a solution u — {ut.x', {t, x) G [0, 1] x M.'^}. 
Then the solution is unique and admits the Wiener chaos expansion ([3]), where 

Inifn{-,t,x)) ^ / fn{tl,Xi, . . . ,tn,Xn)dWt,,xi ■ ■ - dWt^.x^ 
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is the multiple Wiener integral with respect to W, and /„ S T-LV®^ is the 
symmetric function given by (|4]). (See (7.4) of ITT], or (4.4) of [12], or Proposition 
3.2 of [3J). 

Using and the orthogonality of the Wiener chaos spaces, we conclude 
that for any {t,x) G [0, 1] x R'' and (s, y) £ [0, 1] x W^, 

E[ut.j,Us^y\ = E{ut.x)E{Us^y)+^E[In{fn{-,t,x))In{fn{-,S,y))] 

n>l 

= w{t,x)w{s,y) + ^—^an{t,x,s,y), (11) 

r!>l 

where an{t,x, s,y) (n!)^(/„(-, i, x), /„(•, s, y))-^^®" for any n > 1. Note that 

^y)= I / TT?K*i'Sj)(G't:x,Gs:y)p(Rd)»«dsdt, (12) 

J[nA]^- J[oA^ j=i 



cXyi (i, X, s. 



where ViW^) is the the completion of {lyi;j4 e ^^(M'')} with respect to the 
inner product (1^, Ib)v{vl'') = /a /b /(^ ~ y)dydx, 

n 

Gt-^{xi,...Xn) = w{tp{l),Xp(l))Y[Pt,^, + rj^t,U){Xp(j + l) -Xp^j)) 



Gs;yiyi, ■■■yn) 



^(s^(i),ya(i))nP«-u+i)-^= 



(j)(y<T(j+i) -y<TO)), 



the permutations p and a are chosen such that: 

< ip(i) < tp(2} < ■ ■■ < tp(n) and < s^(i) < s„(2) < ■ ■ ■ < 

tp{n+l) = ti Sa(n+l) — S, X pf^^+l) = ^ ^ud ^^(n+l) = ?/• 

Using (flU]) and (|12p . we conclude that a„(t) admits the representation: 



an{t,x,s,y) = n\ e^" ^ i? 

ii,...,2„>l distinct 



N 



(^t-Tii ,...,t-Ti„ Gs-pi-^ ,---,s-pi^;y)v{R'')S>^dj^.^ .^ 



(13) 



The next result gives a probabilistic representation of the inner product 
appearing in (fT3)) . This result is the analogue of ([6]), needed for the treatment 
of the fractional case. 

Lemma 2.2 For any ti, . . . , t„ G [0,t], si, . . . , s„ G [0,s], x £ and y G M'', 

{Gt-ti,...,t~t„;x, Gs-si,...,s-s„;y)'P{'!Si'')'S'" — 



E 



it - t*,Bl, )w{s ~ s*,Bl, ) n f{Bl^ - Bl ) 
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where t* = max{fi, . . . , f„}, s* = max{si, . . . , s„}, and = (i?j)tg[o.i] and 
= (i3|)jg[o_i] are independent d- dimensional Brownian motions, starting 
from X, respectively y. 

Proof: Let B^ and B^ be independent d-dimensional Brownian motions, start- 
ing from 0. Note that 

I := {Gt-ti,...,t-tn;x,Gs-si,...,s-Sn\y)v(9.<i)®^ = / w{t-tp(i),Xp(i))w{s-Scr(l),ya(l)) 

ri, 

n -^(^J ~ yj)PtM)-tpU+i) (^pU) - a;pO+i))Ps<,o)-s„o+i) iVaU) - yau+i))dydx, 

3=1 

where the permutations p and a of {1, . . . , n} are chosen such that 

< tp(„) < . . . < tp(i) <t and < s^(„) < . . . < s^(i) < s, 
tp{n+i) = Sa{n+i) = 0, Xp(^n+i) = X, y^(„+i) = y. We use the change of variables 

^pU) ~ ^pU+1) = ^n+l-j, ycr(j) " VaU+l) = Wn+l-j, j = l,...,n. 

Note that Xp(^j) = x + jyk=i~^ ^k, i-c Xj = x + J2k=i^'' ^k- We get: 
/= / w{t - tpn),x + Yzk)w{s - s„m,y + S2wk) 

n / + 23 zk)-{y+ ^'') 

j=l \ k=l fc=l 

n n 
n^^t.O)-«PO+i) i^n+l-3)J{Ps,u)-^^i3+^) {Wn+l-j)dwdz. 

We now use the fact that pt-s{x)dx is the density of the increment of a d- 
dimcnsional Brownian motion over the interval [s, t], these increments over dis- 
joint intervals are independent, and the Brownian motions B^ , B^ are indepen- 
dent. Therefore, we replace the integral over M^"** by the expectation , 
and the variables Zk,Wk by B} — B} respectively , , — 

, ^ , for all fc = 1, . . . , n. Note that, for any m = 1, 



...,n 



YiBl -Bl )=B} 

tp(r. + l-fc) tp(n + l-fc + l) tp(n + l-m) 

fe=l 

V^/p2 — R2 ) = 

/ So-(n + l-fc) S<^(„+i_fe+i) y S<r(n+l-m)' 

fe=l 
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Hence, 

n 

□ 

Conclusion of the Proof of Theorem [TTT} Since equation ([7]) has a solution, 
this solution is unique and admits the Wiener chaos expansion The result 
follows from ([TTl), ^ and Lemma O □ 
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